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ABSTRACT 

.^L ' Context. Many active galactic nuclei exhibit X-ray features typical of the highly ionized gas called "Warm Absorber" (WA). Such a material 
, appears to be stratified, displaying zones of different density, temperature, and ionization. In this paper, we investigate the possibility of 
^—-j ■ modelling the WA gas in NGC 3783 as a single medium in total pressure equilibrium. 

^) Aims. Our goal is to demonstrate that the WA can be well modelled assuming constant total pressure, in contrast to the current descriptions that 
, are based on the presence of multiple regions, each in constant density. The assumption of total pressure equilibrium yields a more physical 
description of the WA, resulting in the natural stratification of the ionized gas, and providing an explanation for the presence of lines from 
different ionization states, as observed in WA spectra. 

Methods. We have used the photoionization code TITAN, developed by our team, to compute a grid of constant total pressure models with the 
O ' purpose of fitting the WA in NGC 3783. We have compared our models to the 900 ks Chandra spectrum of NGC 3783 and to previous studies 
where the WA was described by multiple zones of constant density. 
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■ Results. In the case of NGC 3783, the WA features can be well reproduced by a clumpy, ionized gas with cosmic abondances, ionization 



C3 



parameter % = 2500 erg cm s , column density Nh = 4 10 22 cirT 2 , and constant total pressure. 

Conclusions. We have shown that the WA in NGC 3783 can be modelled by a single medium in total pressure equilibrium; this is probably 
the case for other WAs currently described by multi-zone, constant density models. In addition, our work demonstrates that the TITAN code is 
well adapted to the study of the WA in active galactic nuclei, opening new prospects for the use of TITAN by a larger community. 
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1. Introduction 

Many Active Galactic Nuclei (AGN) exhibit important X-ray 
absorption features caused by the presence of highly ionized 
gas located in the line-of-sight of the central continuum; this 
material is called "Warm Absorber" (hereafter WA). The first 
observations of WA gas in AGN were reported by Halpern et al. 
(1984) in the Einstein Observatory spectrum of MR 2251-178, 
a quasar displaying a large absorption feature around 1 keV; 
such a feature is consistent with the presence of gas photoion- 
ized by the hard X-rays produced near the AGN central engine. 

Early ASCA observations have revealed the presence of in- 
trinsic absorption in ~50% of type 1 Seyferts. With the ad- 
vent of space X-ray observatories such as Chandra and XMM- 
Newton, which carry aboard high-resolution grating spectro- 
graphs, an important set of high quality data became avail- 
able; evidence for the presence of warm gas was found in 
many other types of AGN, e.g. Seyfert 2s (Kinkhabwala et al. 
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2002), Narrow Line Seyfert Is (e.g. Matsumoto et al. 2004), 
BAL QSOs (Guainazzi et al. 2001 ; Grupe et al. 2003), and even 
some BL Lacs (e.g. Sambruna et al. 1997), although opinions 
diverge for these objects (Blustin 2004). In Seyfert Is, the WA 
spectra revealed the presence of tens of absorption lines, cov- 
ering a wide range of ionization states and blueshifted by a few 
hundreds to thousands of km s _1 (an indication that the absorb- 
ing material is outflowing). In type 2 AGN, the data have shown 
the presence of emission lines. 

Despite the undeniable improvements in our knowledge of 
the WA, some important issues remain a subject of debate, 
namely: (i) the geometry (multi-zone gas? clumps or flow?) 
and location of the WA, (ii) the physical conditions of the ab- 
sorbing/emitting gas, and ( Hi) the importance of the WA in the 
energetics of the AGN. Trying to solve these questions requires 
not only high quality observations, but also an adequate mod- 
elling of the X-ray data through the use of reliable photoioniza- 
tion codes, which calculate the full radiative transfer and allow 
for a good physical description of the ionized gas. 
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The WA gas appears to be stratified, displaying multiple 
zones of different density, temperature, and ionization (e.g. 
Krolik & Kriss 2001 ; Rozariska et al. 2006). Until now, the WA 
in AGN has been modelled assuming the presence of several re- 
gions, each region being modelled in constant density and hav- 
ing a different ionization parameter 1 (^, in units of erg cm s _1 ), 
column density (Nh, in units of cm -2 ), etc. However, a strati- 
fied medium can be obtained naturally, if one assumes the gas 
to be in pressure equilibrium. Such a medium could explain 
the presence of lines from different ionization states and ac- 
count naturally for the other properties of a model composite 
of several constant density regions. 

We have addressed the above-mentioned points through the 
study of the WA in NGC 3783, for which unrivaled quality 
Chandra/HETGS data (Kaspi et al. 2002) are available in the 
archives. There is a short discussion of the WA in NGC 3783 
in the next section. We have modelled the data using the pho- 
toionization code TITAN, briefly described in Sect. 3. Our first 
results are given in Sect. 4, and the conclusions in Sect. 5. 

2. The Warm Absorber in NGC 3783 

NGC 3783 is a bright (V-13.5), nearby (z = 0.0097), 
Seyfert 1.5 galaxy observed in the Optical, UV and X-rays. 
The WA in this object has been discussed by several authors 
(e.g. Kaspi et al. 2001,2002; Netzer et al. 2003; Krongold 
et al. 2003; Behar et al. 2003) based on Chandra and XMM- 
Newton data. These studies seem to agree on the presence of a 
multi-zone gas constituted by three or more regions of differ- 
ent ionizations. In terms of the kinematics of the WA, two or 
more velocity systems have been identified in Chandra obser- 
vations; they are compatible with those observed in UV spec- 
tra. A single outflow velocity system (va ow ~ 600-800 kms -1 ) 
seems enough to describe XMM-Newton observations, which 
have lower spectral resolution. There is no consensus on a pos- 
sible correlation between the velocity shifts, or the FWHMs, 
with the ionization potentials of the ions. 

Although the WA in NGC 3783 has been the object of 
many studies, these have assumed either a constant density (e.g. 
Netzer et al. 2003), or a dynamical state (Chelouche & Netzer 
2005) for the modelling. In addition, they all require multiple 
zones of different density, temperature, and ionization, which 
are invoked to explain the large span in ionization observed in 
the WA spectrum. Furthermore, when plotted on the S-curve of 
thermal equilibrium log (T) vs. log (E,/T) (where T is the tem- 
perature of the medium and § is the ionization parameter), these 
zones lie on a vertical line of roughly the same gas pressure, 
suggesting that a single medium in pressure equilibrium could 
account for the WA in this object. A physical modelling of the 
gas through adequate photoionization codes remains necessary 
to confirm or to invalidate such an assumption, that being the 
main goal of this paper. 

We note that, for relatively thin slabs, the radiation and gas 
pressure remain constant, the temperature does not vary much 

1 The ionization parameter is defined as L/n^R 2 , where L is the 
source's bolometric luminosity, « H is the hydrogen number density at 
the illuminated face of the medium, and R is the distance from the WA 
to the illuminating source. 



within the slab, and constant pressure models can thus be ap- 
proximated by a constant density model. However, for thicker 
media, the ratio between the gas and the radiation pressure 
varies across the slab, inducing variations in the density and 
the temperature. In such cases, it is important to take the total 
(gas+radiation) pressure into account. 

3. Modelling the WA with the TITAN code 

We have used the photoionization code TITAN, developed by 
our team (Dumont et al. 2000,2002; Collin et al. 2004), to 
model the 900 ks C/zanc/ra/HETGS observations and to con- 
strain the physical conditions of the WA gas in NGC 3783. 

TITAN is well-suited both for the study of optically thick 
(Thomson thickness up to several tens) and thin media, such as 
the WA. It computes the gas structure in thermal and ioniza- 
tion equilibria, giving as output the ionization and temperature 
structures, as well as the reflected, emitted outward, and ab- 
sorbed spectra. One advantage of TITAN is that it treats the 
transfer of both the lines and the continuum using the ALI 
(Accelerated Lambda Iteration) method, which very precisely 
computes line and continuum fluxes. Our line transfer treat- 
ment leads to a longer computation time than with other ap- 
proximate methods; we thus use fewer lines than those avail- 
able in other codes. Our atomic data include ~ 10 3 lines from 
ions and atoms of H, He, C, N, O, Ne, Mg, Si, S, and Fe; they 
omit features like the Unresolved Transition Array (UTA) and 
the inner shell transitions (except the iron K lines). In the case 
where the radiation pressure and the losses from the lines are 
only a few tens percent of those from the continuum, as is the 
case in this study, the missing lines do not change the ionization 
and thermal structure much (Chevallier et al. 2006). 

The total (gas + radiation) pressure is computed after the 
plasma temperature has been determined by the thermal bal- 
ance equation. The pressure is kept fixed within the medium, 
thus determining the density profile. The model is iterated un- 
til convergence (see Dumont et al. 2000, 2002). In the presence 
of an external, strong, irradiating source, such as the hard X- 
rays ionizing the WA material, the irradiated medium develops 
an ionization instability (Krolik et al. 1981). If the medium is 
in pressure equilibrium, it leads to a temperature discontinuity 
(e.g. Ko & Kallman 1994) with at least three states of ther- 
mal equilibria (high-, intermediate-, and low-temperature), the 
middle one being unstable. Our approach to the thermal insta- 
bility problem consists of adopting an equilibrium temperature 
value between the "hot" and the "cold" stable solutions. Such 
a method provides a unique solution (for a more complete dis- 
cussion of the thermal instabilities with TITAN, see Sect. 4.1.3 
in Rozanska et al. 2002). 

Our models assume a ID plane-parallel geometry, with the 
WA medium being illuminated on one side by an incident con- 
tinuum. We assume the incident radiation field to be concen- 
trated in a very small pencil-like shape centered on the nor- 
mal direction. The models are parameterized by the ioniza- 
tion parameter the column density Nh, the hydrogen number 
density n H , and the incident flux Spectral Energy Distribution 
(SED). We can also consider a non-zero turbulent velocity v tur b- 



A.C. Goncalves et al.: A new model for the Warm Absorber in NGC 3783: a single medium in total pressure equilibrium 



3 



Concerning NGC 3783 specifically, we note first that the 
turbulent velocity deduced from the observations (e.g. Kaspi et 
al. 2000) is highly supersonic compared to the thermal veloc- 
ity (v t h) in the absorbing medium. It is possible to address this 
situation either (z) by taking into account the strong dissipa- 
tion occuring through shock waves in the supersonic medium 
(this leads to a higher temperature than in a purely photoion- 
ized case), or (ii) by assuming a photoionized, inhomogeneous 
medium with several clumps lying in the line-of-sight. These 
gas clumps have an average relative velocity, that we identify 
as Vturb- We thus assume the WA in NGC 3783 to consist of 
such gas blobs embedded in a surrounding, hotter medium. The 
clumpy WA medium is assumed to be in total pressure equilib- 
rium. We include the gas and the radiation pressure in the total 
pressure, but not the turbulent pressure resulting from taking 
a non-zero v tur b, since it is a quantity exterior to the clumps. 
One should, however, take the effects of v tul -b into account when 
computing the radiative transfer; this has been done for all our 
models. 

The turbulent velocity acting on the line transfer is typically 
one order of magnitude larger than the thermal velocity of the 
heavy elements. The optical depth being inversely proportional 
to v tul b, this difference in velocity can be accounted for if each 
clump has a column density < Nn/n (where A^h corresponds to 
the column density of the clumpy system and n is the number of 
clumps in the line-of-sight). From the observer's point of view, 
and for the purposes of our modelling, these clumps can be 
treated just like a continuous medium with column density A^h, 
as each clump receives the radiation transmitted - and returned 
- by the other ones, while the hot surrounding medium does not 
have any effect on the radiation transfer. 

We have computed a grid of 16 constant total pressure mod- 
els to fit the WA in NGC 3783, covering the combinations of 
the 4 possible values of the ionization parameter (2000 < 'E, < 
3500) and the total column density (3 1 22 < N H < 6 10 22 ). 
The «h was set to 10 5 cm 4 (Netzer et al. 2003) and the v tur b to 
150 kms _1 (Kaspi et al. 2000). The incident SED used, given 
in Kaspi et al. (2001), covers the 0.2eV-400 keV range. The 
models were computed using two different sets of abundances 
(Allen 1973; Netzer et al. 2003). 

4. Results and discussion 

Our study shows that the WA can be modelled under total 
pressure equilibrium conditions, and provides a best model 
with '£, = 2500 erg cm s _1 , Nu = 4 10 22 cm -2 , and cosmic (Allen 
1973) abundances. While the results are stable for a «h up to 
10 12 cirT 3 , the turbulent velocity of the clumpy medium seems 
to play a non-negligible role in the modelling. Our best model 
gives a reasonably good fit to the observed data, both for the 
continuum (following its overall shape up to 10 keV and re- 
producing the O vii and O vm edges) and the lines (both from 
high and low ionization). The observed and modelled spectra, 
as well as the corresponding EWs, will be presented in a forth- 
coming paper (Goncalves et al., in preparation), in which the 
effects of using different abundances and v tur b values will also 
be discussed. 
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Fig. 1. The solid line represents the temperature profile result- 
ing from our modelling of the WA as a single medium in total 
pressure equilibrium. In this case, the temperature discontinu- 
ities arise naturally. The dashed line corresponds to a single- 
zone WA with the same parameters, only modelled under con- 
stant density. 

For all the models in our grid, irradiation by the incident 
SED resulted in the formation of a hot surface zone, with an 
almost constant temperature, followed by a rapid drop by one 
order of magnitude at a certain depth, and a second temperature 
drop deeper in the slab. In Fig.^ we gi ye the temperature pro- 
files for a constant density and a constant total pressure model, 
computed assuming the same parameters. We observe that, in 
the case of a constant density model, the temperature remains 
more or less constant along the medium. Thus, if a composite 
(i.e. multi-component) constant density model is assumed, as 
in Netzer et al. (2003), a relatively flat temperature profile is 
found for each component, with a different temperature value 
for each one. When a single medium in constant total pressure 
is assumed, the temperature drops arise naturally; such varia- 
tions in temperature are accompanied by ionization fronts, re- 
sulting in different ionization fractions for the three sections 
(Tl, T2, and T3) visible in the plot. 

Our results can be compared to those of Netzer et al. (2003), 
who provide a detailed modelling of the WA in NGC 3783; 
these authors found three constant density components with 
different ^ and A^h values. A full comparison between the multi- 
component model of Netzer et al. (as computed by TITAN) 
and our constant total pressure model will also be presented 
in a forthcoming paper. Figure |2]displays an example with the 
computed ionic column densities deduced from our best model 
and from the composite model of Netzer et al. vs. the "ob- 
served" ones determined through the curve of growth analysis 
in Netzer et al. (2003). The figure illustrates the good agree- 
ment between a single medium in constant total pressure and 
the three-zone model in constant density showing that, within 
the uncertainties, both models should lead to similar absorption 
spectra; this comparison used the same incident SED and abun- 
dances. Although Fig.|5]only shows 15 ions (for comparision 
with Netzer et al.), our model provides ionic column densitites 
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Fig. 2. Computed ionic col- 
umn densities (in cm -2 ) for the 
Netzer et al. (2003) compos- 
ite constant density model (left 
panel) and for our constant total 
pressure model (right panel) vs. 
the "observed" column densities 
determined through a curve of 
growth analysis. The observed 
O vi column density value cor- 
responds to a lower limit, only. 
It was therefore replaced by the 
Netzer et al. computed value, 
which explains the absence of 
an error bar for this ion. The 
dotted line shows the diagonal 
on which the calculated values 
equal the observations. 



for all the other ions in the atomic data, which include the lower 
ionization species that give rise to absorption lines in the UV. 

Based on our best model results, and on the object's 
bolometric luminosity (L~2 10 44 erg s _I ) and black hole mass 
(Mbh~3 10 7 M q ) (Peterson et al. 2004), we have calculated 
some quantities related to the WA. Assuming «h = 10 5 cm -3 , 
the sum of the thickness of all the clumps in the line-of- 
sight amounts to AR c i ~2 10 17 cm (a higher «h would imply 
a smaller AR C \). Knowing AR C \, and to keep M out /MEdd ^ 1 
(where M out and Medd are the outflow and Eddington mass ac- 
cretion rates), the volumic factor f wo \ should be < 0.5 (see Eq. 4 
in Chevallier et al. 2006). As f m \ corresponds to the size ratio 
of the "clumpy" to "whole WA" media (A/? c i/A/?waX we can 
estimate A/?wa as being >4 10 17 cm. Based on the previously 
given definition of we calculate the distance to the WA to be 
R < 10 18 cm (for n H = 10 s cirT 3 ). 

Given these values, the replenishment timescale (/?/vfl ow ) 
in NGC 3783 is of the order of 5 10 2 yrs, and the dynamical 
timescale, several times higher. A satisfactory WA geometry 
could be obtained with f m \ ~ 0.2 and n > VturbMh ~10. The geo- 
metrical thickness of a single clump should thus be < 2 10 16 cm, 
resulting in a sound crossing time of only ~ 50 yrs; if n > 10, 
this time will be reduced in proportion. In such a picture, the 
replenishment and dynamical times are larger than the sound 
crossing time in each clump. In consequence, the clumps can 
adjust promptly to any change in the pressure equilibrium of 
the whole system. We note that these timescales are still much 
longer than the variation time of the illuminating flux; it could 
thus be objected that the gas in NGC 3783 cannot reach a state 
where it is in total pressure equilibrium with the illuminating 
source. We suggest that, unless the flux variations are very 
large, the medium would adopt a "quasi-pressure equilibrium" 
corresponding to an incident flux averaged over a long time. 
The flux and spectral variations should thus induce changes 
in the temperature and ionization equilibrium, but the density 



structure would stay about the same. The way the spectrum 
would be modified will be discussed in a forthcoming paper. 

5. Conclusions 

1 . We have shown that the WA in NGC 3783 can be modelled 
by a single medium in total pressure equilibrium. 

2. Our grid of models has provided a best result for § = 2500 
erg cm s _1 , Nh =4 10 22 cm -2 , and cosmic abundances. 

3. Our work opens new prospects for the future use of the 
TITAN code by the community through larger grids of con- 
stant total pressure models to be made available via XSPEC 
and/or via Virtual Observatory facilities. 
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